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FLUID PROPERTIES MEASUREMENTS 
USING WAVELENGTH MODULATION 
SPECTROSCOPY WITH FIRST HARMONIC 
DETECTION 

CROSS-REFERENCE TO RELATED 
APPLICATIONS 

This application claims priority to and the benefit of the 
filing of U.S. Provisional Patent Application Ser. No. 61/636, 
937, entitled “Fluid Properties Measurements Using Wave- 
length Modulation Spectroscopy with First Harmonic Detec- 
tion”, filed on Apr. 23, 2012, and the specification and claims 
thereof are incorporated herein by reference. 

STATEMENT REGARDING FEDERALLY 
SPONSORED RESEARCH OR DEVELOPMENT 

This invention was made with government support under 
Contract No. FA8650-11-M-2146 awarded by the U.S. Air 
Force and No. NNX13CS03P awarded by the National Aero- 
nautics and Space Administration (NASA). The government 
has certain rights in the invention. 

INCORPORATION BY REFERENCE OF 
MATERIAL SUBMITTED ON A COMPACT DISC 

Not Applicable. 

COPYRIGHTED MATERIAL 

Not Applicable. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention (Technical Field) 

The present invention relates to methods and apparatuses 
for determining the velocity, temperature and pressure of 
gases in supersonic inlets of scramjet engines, and any other 
gas flows requiring the knowledge of pressure, velocity and 
temperature. 

2. Description of Related Art 

Note that the following discussion refers to a number of 
publications by author(s) and year of publication, and that due 
to recent publication dates certain publications are not to be 
considered as prior art vis-a-vis the present invention. Dis- 
cussion of such publications herein is given for more com- 
plete background and is not to be construed as an admission 
that such publications are prior art for patentability determi- 
nation purposes. 

Turbine-based combine cycle (TBCC) and rocket-based 
combined cycle (RBCC) are airbreathing propulsion tech- 
nologies for two-stage-to-orbit vehicles that operate at costs 
lower than conventional methods based solely on rockets and 
scramjets. TBCC is applicable in the flow regime from sub- 
sonic to supersonic and not suitable for orbital insertion, but 
is the choice for a high speed passenger transport (Sippel, M. 
and Okai, K., “Preliminary Definition of a TBCC Propulsion 
System for a Mach 4.5 Supersonic Cruise Airliner,” ISABE 
2007-1204, International Society for Air Breathing Engines). 
RBCC is suitable for hypersonic to rarefied flow regimes. 
RBCC based vehicles are capable of reaching high altitudes 
and space with pure rocket operation (Kloesel, K. J., Ratnay- 
ake, N. A., and Clark, C. M., “A Technology Pathway for 
Airbreathing, Combined-Cycle, Horizontal Space Launch 
Through SR-71 Based Trajectory Modeling,” AIAA 2011- 
2229, 17th AIAA International Space Planes and Hypersonic 
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Systems and Technologies Conference, San Francisco, 
Calif.). The need for the non-intrusive measurement of veloc- 
ity, pressure and temperature profiles at the inlet and exhaust 
flows is obvious when supersonic flow speeds are present. 
5 Traditional probes used for temperature, pressure and veloc- 
ity disturb the flow sufficiently that corrections, based on a 
priori measurements and models, are needed to ascertain the 
local flow properties. Moreover, wall-based measurements of 
pressure do not provide information on local nonuniformities . 
111 Even in uniform flows, pressure in the core flow is determined 
via corrections using models with the knowledge of the adja- 
cent wall pressure. Temperature and velocity measurements 
based on optical methods such as particle imaging velocim- 
15 etry (PIV) and planar laser-induced fluorescence (PLIF) are 
excellent fora small region of interest but not easily amenable 
to a full profile of the inlet and exhaust flows. Accurate infor- 
mation of inlet flow properties is important for gauging the 
mass capture and properly tuning the fuel flow for efficient 
20 combustion. Flow properties of exhaust flows are important 
for assessing the thrust generated by the combustion system 
(Drummond, J. P., Cockrell, C. E., Jr., Pellett, G. L., et al., 
“Hypersonic Airbreathing Propulsion — An Aerodynamics, 
Aerothermodynamics, and Acoustics Competency White 
25 Paper,” NASA/TM-2002-21 1951). 

Velocity and temperature measurements via tunable diode 
laser absorption spectroscopy (TDLAS) are well established 
for measuring gas temperatures and velocities in combustor 
flows. Both direct absorption and wavelength modulation 
30 spectroscopy (WMS) methods exist to make these measure- 
ments (Hanson, R. K., “Applications of quantitative laser 
sensors to kinetics, propulsion and practical energy systems,” 
Proc. Combust. Inst. 33, 1-40 (2011); Li, F., et al., “Uncer- 
tainty in velocity measurement based on diode-laser absorp- 
35 tion in nonuniform flows,” Appl. Opt. 51, 4788-4797 (2012); 
Barhorst, T., Williams, S., Chen, S.-J., Paige, M. E., and 
Silver, J. A., “Development of an In Flight Non-intrusive 
Mass Capture System,” AIAA Paper No. 2009-5067, 45th 
AIAA/ASME/SAE/ASEE Joint Propulsion Conference & 
40 Exhibit (2009); Silver, J. A., “Frequency Modulation Spec- 
troscopy for Trace Species Detection Theory and Compari- 
son Among Experimental Methods,” Appl. Opt. 31, 707-717 
(1992); Bomse, D. S., Stanton, A. C., and Silver, J. A., “Fre- 
quency Modulation and Wavelength Modulation Spec- 
45 troscopies: Comparison of Experimental Methods Using a 
Lead-Salt Diode Laser,” Appl. Opt. 31, 718 (1992)). These 
studies always require the recording of the full absorption 
spectrum, and for temperature measurements often require 
two lasers. For high speed (>1 kHz) transient measurements, 
50 the electronic requirements for high precision and sensitivity 
become cumbersome at best, where data acquisition rates and 
demodulation electronics pushbeyond MHz frequencies. The 
analyses of these measurements can be quite complex, requir- 
ing measurement of a host of spectroscopic parameters and 
55 their dependence on temperature and pressure in order to 
retrieve the local gas temperature. This makes the sensor 
systems inflexible and suitable for only one gas, whereas 
multiple gases may be of interest. All of the aforementioned 
measurements apply to velocity or temperature. There have 
60 been no reported methods to optically measure local pressure. 

Hanson (U.S. Pat. No. 5,178,002) discussed a system and 
method for determining the value of thrust non-intrusively 
from products of combustion of a jet engine. Both line-of- 
sight absorption spectroscopy and laser-induced fluorescence 
65 methods are discussed for measuring pressure, velocity, and 
temperature of a gas. Full absorption spectra are used in the 
measurements. 
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Williams, S., et al., “Diode Laser Diagnostics of High 
Speed Flows,” AFRL-PR-WP-TP-2007-204, Air Force 
Research Laboratory, Propulsion Directorate, October 2006, 
discussed the use of direct absorption spectroscopy with full 
line shapes and multiplexing of multiple lasers for the mea- 
surement of mass capture in hypersonic inlets and isolators. 

Lindstrom, C., et al., “Diode Laser Absorption Tomogra- 
phy of 2D Supersonic, Flow,” 43 rrf AIAA/ASME/SAE/ASEE 
Joint Propulsion Conference & Exhibit, 8-1 1 Jul. 2007, Cin- 
cinnati, Ohio, examined the combined use of tomography 
reconstruction techniques and direct absorption spectroscopy 
with full line shape to probe supersonic flows for determining 
flow properties and species concentrations. 

Chen, S.-J., et al., “Laser-Based Mass Flow Rate Sensor 
Onboard HIFiRE Flight 1,” 45* AIAA/ASME/SAE/ASEE 
Joint Propulsion Conference & Exhibit, 2-5 Aug. 2009, Den- 
ver, Colo., used wavelength modulation spectroscopy with 
second harmonic detection and full line shape to measure gas 
velocity in high speed flows. The measurements were carried 
onboard a sounding rocket reaching a speed of Mach 8. 

Wheatley, B., “Tunable Diode Laser Absorption Spectros- 
copy of Hypersonic Flows,” 28* International Congress of 
the Aeronautical Sciences,” ICAS, 2012, reviewed the 
method of tunable diode laser absorption spectroscopy for the 
measurement of flow properties and species concentrations in 
hypersonic and scramjet flows. Direct absorption spectros- 
copy and wavelength modulation spectroscopy using first and 
second harmonic detection, using full line shape, were com- 
pared. 

Girouard, R., et al. (U.S. Pat. No. 8,265,851) disclosed a 
method for controlling the performance of engines via the 
measurement of a wavelength-dependent parameter. Based 
on the measured parameter, a needed adjustment to the com- 
bustion event is determined and physically applied to the 
engine to enhance engine performance. Direct absorption 
spectroscopy and wavelength modulation spectroscopy were 
identified as some of the measurement methods to obtain 
velocity, temperature, and species concentrations, all using 
full line shapes. 

The present invention improves on the art by, e.g., using an 
optical velocity measuring system that comprises a light 
source from a diode laser and two detectors to monitor two 
crossed-beams paths. Alternatively, pressure can be mea- 
sured using a single path across the flow. The present inven- 
tion provides the following advantages over the state-of-the- 
art sensors: (1) not requiring the need to scan the full 
absorption line shape; (2) no need to fit the data to obtain peak 
positions of measured absorption line; (3) no modeling of line 
shapes to account for residual amplitude modulation; (4) 
capable of achieving high-bandwidth measurements of 10 
kHz or higher; (5) measuring only two point values per 
probed path; (6) inherently self-calibrating; (7) no contribu- 
tions from interference fringes (etalons) which show up only 
as dc offset; and (8) high-sensitivity using wavelength modu- 
lation spectroscopy with If detection. 

BRIEF SUMMARY OF THE INVENTION 

The present invention is of a method (and concomitant 
apparatus) for non-intrusively obtaining the local velocity, 
pressure, and temperature of a gas, comprising: emitting two 
beams of laser light directed at two different angles across a 
flow of the gas; separately detecting the two beams; deter- 
mining first harmonic line shapes from the detected beams 
using wavelength modulation spectroscopy and a pair of 
single point measurements; calculating gas velocity and pres- 
sure from Doppler shifts of the line shapes; and calculating 
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gas temperature from slopes of the line shapes; wherein 
wavelength of the two beams is alternated between two fixed 
wavelengths; and using a reference cell to provide zero-ve- 
locity reference and pressure reference. In the preferred 
5 embodiment, a square wave is employed, most preferably 
wherein frequency of the square wave is at or greater than 
desired detection bandwidth. Scanning full absorption line 
shapes is not required. Gas velocity is determined indepen- 
dent of gas density changes. Calculating gas temperature 
10 comprises employing a mole fraction of the gas and/or mea- 
suring second harmonic peak signals. The invention prefer- 
ably operates at a bandwidth of up to about 40 kHz, obtains 
velocity at a resolution of about 10 m/s fluctuation or better, 
and obtains pressure at about plus or minus 12 torr at 2 
15 atmospheres. 

Further scope of applicability of the present invention will 
be set forth in part in the detailed description to follow, taken 
in conjunction with the accompanying drawings, and in part 
will become apparent to those skilled in the art upon exami- 
20 nation of the following, or may be learned by practice of the 
invention. The objects and advantages of the invention may be 
realized and attained by means of the instrumentalities and 
combinations particularly pointed out in the appended 
claims. 

25 

BRIEF DESCRIPTION OF THE SEVERAL 
VIEWS OF THE DRAWINGS 

The accompanying drawings, which are incorporated into 
30 and form a part of the specification, illustrate one or more 
embodiments of the present invention and, together with the 
description, serve to explain the principles of the invention. 
The drawings are only for the purpose of illustrating one or 
more preferred embodiments of the invention and are not to 
35 be construed as limiting the invention. In the drawings: 

FIG. 1 is a schematic diagram illustrating the apparatus for 
measuring pressure and/or temperature using a single optical 
path that is perpendicular to the flow path; 

FIG. 2 is a schematic diagram illustrating the apparatus for 
40 measuring velocity, pressure, and temperature using two opti- 
cal paths angled in opposite direction to the flow path; 

FIG. 3 is a schematic diagram illustrating the velocity- 
induced and pressure-induced shift of wavelength with 
respect to wavelength at zero velocity and at reference pres- 
45 sure using wavelength modulation spectroscopy with If 
detection; 

FIG. 4 is a schematic diagram illustrating the method of 
Doppler shift for measuring velocity using wavelength modu- 
lation spectroscopy with 1 f detection and two points along the 
50 linear portion of the curve; and 

FIG. 5 is a schematic diagram illustrating the method for 
measuring pressure using wavelength modulation spectros- 
copy with If detection using two points along the linear 
portion of the curve. 

55 

DETAILED DESCRIPTION OF THE INVENTION 

The present invention is of a method of and apparatus for 
monitoring non-intrusively gas velocity, pressure and tem- 
60 perature in combustion systems and flow devices. It is capable 
of measuring these parameters (velocity, pressure, and tem- 
perature) without the need to record the full absorption spec- 
tra of the probed chemical species. It is capable of monitoring 
these parameters at faster speeds compared to traditional 
65 methods which record the full absorption spectra. It is capable 
of measuring these parameters simultaneously using the same 
optical paths and setups. It is capable of providing a means to 
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communicate with the systems or devices for the purpose of 
active-control to enhance efficiencies. It is capable of inter- 
rogating multiple locations simultaneously within the com- 
bustion systems and flow devices. 

In the preferred apparatus for pressure and temperature 5 
measurements, as illustrated in FIG. 1 , a computer 4 (or other 
waveform generating source) generates a slow square wave 
(i.e., low frequency) and fast sine wave (i.e., high frequency). 

A laser controller 6 uses the slow square wave and fast sine 
wave to control a laser source 8 . A beam splitter 10 directs the to 
laser beam into two paths. One path sends the light into a 
line-locking cell 1 6 and another portion of the light into a light 
collimator 12 for focusing. Light exits the beam collimator 12 
and traverses the probed region 2 and is collected by a detec- 
tor 14 . Light traversing the line-locking cell is collected by 15 
another detector 18 . Hie AC-component of signal from the 
detector 18 is processed by a phase-sensitive detector 20 and 
sent to the computer 4 for additional processing. The AC- 
component of light collected by detector 14 is also processed 
by a phase-sensitive detector 22 . The DC-component of the 20 
signal from detector 14 is processed directly by the computer 
4 . 

In the preferred apparatus for velocity measurements, as 
illustrated in FIG. 2 , a computer 4 (or other waveform gen- 
erating source) generates a slow square wave (i.e., low fre- 25 
quency) and fast sine wave (i.e., high frequency). A laser 
controller 6 uses the slow square wave and fast sine wave to 
control a laser source 8 . A beam splitter 10 directs the laser 
beam into three separate paths. Two light collimators 12 and 
14 are used to direct the laser beam into two measurement 30 
paths that are crossed by an angle 20 to probe air flow 2 . Two 
detectors 16 and 18 are used to collect signals from these two 
measurement paths. Phase-sensitive electronics 26 and 28 are 
used to collect the measurement AC-signals. The DC-signals 
are processed directly by a computer 4 . A third beam is 35 
directed through a line-locking cell 20 and onto detector 22. 
The AC-signal from detector 22 originating from line-locking 
cell 20 is collected by phase-sensitive detection electronics 24 
and processed by computer 4 . 

In the preferred embodiment the beam injection system 40 
comprises a focusing lens for collimating the laser output 
beam from preferably a fiber optic cable, preferably matched 
for the operating wavelength. The detection system com- 
prises a photodetector, preferably matched for the operating 
wavelength. For pressure measurements, the coupler 10 in 45 
FIG. 1 comprises a fiber splitter having 90 percent and 10 
percent output legs. For velocity measurements (FIG. 3), the 
first coupler comprises a fiberized splitter having 45, 45, and 
1 0 percent output legs. Alternatively, this 3-way splitter could 
be replaced with two splitters, the first having 90 percent and 50 
10 percent output legs. The 90 percent leg is subsequently 
split by a second 50-50 splitter to provide intensity to the 
crossed beams in the flow. Any combination of splitting laser 
power using 2-way and 3-way splitters can be used. 

Advantages of this approach include measurement band- 55 
width of 1 0 kFlz or faster, high-sensitivity measurements with 
absorbance of 1 part 1 ,000 at 1 0 kFlz are possible, inherently 
self-calibrating step, and insensitivity to optical interferences 
like etalons. The method uses a single point measurement 
without the need to scan the full absorption line shape as it is 60 
typically done in TDLAS based measurements. 

Traditionally, high sensitivity detection using optical 
absorption uses wavelength modulation spectroscopy. The 
main features of WMS are that spectral noise is reduced by 
performing signal detection at frequencies, high enough to 65 
avoid laser excess (1/f) noise fluctuations in the laser output 
power (>10 kFlz) and that the measurements are made against 


6 

a zero-valued baseline. To perform WMS, a small high-fre- 
quency sinusoidal modulation is superimposed on the diode 
laser injection current. This current modulation produces a 
modulation of the laser wavelength at the same frequency 
since the laser wavelength is tuned by changing the current. 
Absorption by the target gas converts the wavelength modu- 
lation to an amplitude modulation of the laser intensity inci- 
dent on the detector, adding AC components to the detector 
photocurrent. Phase-sensitive detection of the detector pho- 
tocurrent at a multiple n of the modulation frequency of 
selectively amplifies the AC components and shifts the mea- 
surement from near DC to higher frequencies where laser 
noise is reduced. 

One can also detect at the primary modulation frequency, 
1 f. In this case the WMS line shape looks like the first deriva- 
tive of the direct absorption line shape, where the signal 
crosses the baseline at the absorption peak wavelength. Flere 
the slope efficiency of the laser appears as a dc offset. This 
offset can also be used as a monitor of I 0 . Since the If- WMS 
signal crosses its baseline when the laser wavelength is at line 
center, this signal is often used in a line-locking feedback loop 
to keep the laser scan centered on the peak. A key feature of 
the 1 f line shape is that the line shape near the zero crossing is 
a linear function (i.e., essentially a straight line with no cur- 
vature). 

Typically, one scans the laser wavelength (laser injection 
current) across the spectral line of interest and records the 
absorption spectrum. For very high bandwidth measure- 
ments, this approach is problematic. A typical spectrum 
requires 100-1000 points. To acquire a full spectrum for a 
desired 1 0 kHz bandwidth, one must scan at a minimum rate 
of 10 kHz, with each data point acquired at 1-10 MHz. To 
avoid data aliasing, the modulation frequency must be at least 
three times higher than this rate. A hill spectral scan thus 
requires high speed electronics, and scanning the laser this 
fast reduces its tunability and introduces noise. Data process- 
ing requirements are also more complex and expensive. A 
feature of this invention is that there occurs no scanning 
across the line shape, but the laser wavelength is simply 
alternated between two very closely-spaced wavelengths on 
either side of line center. This feature minimizes the effect of 
etalons, no line shape fitting is required, and high bandwidths 
are readily achieved. 

Gas velocity measurements using absorption employ the 
Doppler effect to relate the frequency shift in the line center 
absorption wavelength to the speed of the gas, as illustrated in 
FIG. 3. A laser beam traverses the gas flow at an angle other 
than 0 degrees, and the observed shift in measured line center 
frequency in wavenumbers v 0 is determined by the compo- 
nent of the velocity vector along the flow direction. The 
shifted line frequency Av c is proportional to the incidence 
angle 0 by, 

Av D =v o (l+Fsin0/c)(l+i iS (P/P o ), (1) 

where V is the mean gas velocity, b s is the shift coefficient 
defined for pressure P 0 , P is the local pressure, c the speed of 
light, and v 0 the peak absorption frequency in wavenumbers. 
This equation accounts for the Doppler shift in a gas flow as 
well as for the effect of pressure-shifting, where the absolute 
line center wavenumber of an absorption peak varies linearly 
with local pressure. Pressure shift coefficients (in units of 
cnrVatm) are tabulated in the HITRAN and GEISA data- 
bases for many common gases (Rothman. L. S., et al., “The 
HITRAN 2008 molecular spectroscopic database,” J. Quant. 
Spectrosc. Radiat. Trans. 110, 533-572 (2009); Jacquinet- 
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Husson, N., et al., “Hie 2009 edition of the GEISA spectro- 
scopic database,” J. Quant. Specrosc. Radiat. Trans. 112, 
2395-2445 (2009)). 

Using a pair of beams crossed at angle 20, as illustrated in 
FIG. 2, both beams shift in the same direction for changes in 5 
pressure, but the Doppler shifts of the crossed beams in a 
moving gas are in opposite directions. For two beams crossed 
at 60 degrees (0=±3O°, the relative separation in line centers 
due to velocity only is 

to 

A v fl ^=v 0 *T/c (2) 

FIG. 3 shows lf-WMS line shapes for such a system of 
crossed beams, referenced to the line shape of a static, fixed- 
pressure gas. Note the reference line shape (dashed line) can 
have a different width (slope) than the sample beams (solid 15 
lines). This figure illustrates how the crossed beams separate 
in wavelength due to their relative (+) and (-) velocity shifts 
in the flow, while both beams move in the same direction if 
their pressure is different than in the reference leg. 

Measuring gas velocity at high bandwidths via the Doppler 20 
shift using a full spectrum is limited by noise, electronics and 
processing constraints. The present innovation eliminates all 
of these complications and permits easy and precise high 
bandwidth measurements. This is illustrated in FIG. 4a where 
the lf-WMS signal is shown for a system that uses (for this 25 
illustration) a 5 kHz detection bandwidth. The lf-WMS sig- 
nal near the center of the line shape varies linearly with 
wavenumber. This effect is explained by the fact that the peak 
of a Voigt line shape can be accurately represented by a 
parabolic function (i.e., y=mx 2 +b). The first derivative (WMS 30 
shape) of this function is a straight line with y' 2mx, where 
x=wavenumber and 2m is the slope. The slope of the 1 f shape 
is related to the gas density. 

Rather than scan the laser across the full absorption spec- 
trum, one could fix the laser at a single wavelength on the if 35 
line shape. As the velocity changes, the line shape would shift 
left or right (as in FIG. 4b) and this signal voltage would 
change in proportion to the velocity shift. As this is now only 
a one-point measurement, velocity can be rapidly deter- 
mined. A full spectrum of many points is not required. How- 40 
ever, as temperature, pressure, or concentration change, so 
does the slope. 

At zero flow velocity, both crossed laser beams produce 
identical line shapes; there is no Doppler shift (FIG. 4a). 
When the gas velocity is non-zero, the two line shapes sepa- 4 - 
rate and shift in opposite directions. Rather than scan across 
these shapes and record full spectra, the approach here simply 
alternates the laser wavelength between two fixed wave- 
lengths on either side of line center and monitors the 1 f signal 
voltages at these two wavelengths. For both crossed beam 50 
channels which are identical here, one measures line shape 
voltages of a=b for the positive square wave step and a'=b' for 
the negative step of the square wave. The difference in signal 
voltage AS between these two values (top and bottom of the 
measured square wave as a function of time) provides the 55 
instantaneous slope, since we know the wavenumber differ- 
ence Av m between the alternating steps. This slope in terms of 
voltage a, a', b, and b' is 

Slope 1 y^(A5'/Av 77l )=(a-a , )/Av 77I =(6-6 , )/Av ra . Units of 60 

volts* cm- 1 (3) 

Now as the velocity becomes non-zero (FIG. 4b), the 
observed square wave voltages for the two beams separate 
vertically, providing a measure of the velocity shift Av n (FIG. 
4b), where 65 

( 4 ) 
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Combining Eqns. (2)-(4) with a previously measured cali- 
bration of Av m , one finds the mean gas velocity is 

F=A v„*c*(a-i)/(v 0 *(a-a')) (5) 

If, for fixed V, the slope of the line shape were to vary due 
to a change in pressure, temperature, or concentration of the 
gas, the individual square wave voltage ranges (a-a' or b-b') 
would change, but the (+) and (-) velocity beam waveform 
separations (a-b or a'-b') would also scale accordingly, so that 
the resultant velocity would be the same. Therefore, the 
inventive approach simultaneously measures both the slope 
and signal level to measure velocity independent of density 
changes. 

In practice, this measurement is accomplished by imposing 
a small square wave step Av„, on the laser current so that the 
measurement wavelength alternates above and below line 
center (noted by orange vertical lines). The frequency of this 
alternation should be at or greater than the desired detection 
bandwidth, so as to capture all fluctuations in density that 
might be correlated to velocity shifts. For example, if one 
desired a detection bandwidth of 5 kHz, the step frequency 
could be 10 kHz, and the analog-to-digital (ADC) acquisition 
rate could be 20 kHz or greater. To be sure that the If modu- 
lation is properly acquired and demodulated, one can set 
f=100 kHz. A helpful feature of this approach is only low 
frequency rf electronics are required, helping assure that 
noise associated with high bandwidth circuits are avoided. In 
addition to these lf-WMS AC-signals from each leg, the 
corresponding DC-signals are also acquired by separate 
ADCs for laser intensity normalization. A significant advan- 
tage of this method is that if the density (P, T or mole fraction) 
changes the slope of the line shape, this effect is immediately 
accounted for in the analysis. Detection bandwidths up to 40 
kHz should be readily achieved. 

The exact location of the absorption peak at zero velocity 
can be found by locking the laser wavelength to a 1 f signal of 
a gas sample in a fixed pressure, static reference cell, and 
using the zero crossing as a feedback signal. This provides a 
reference for zero velocity and for our calibration pressure, 
Po- 

By comparing the square wave signal of one of the velocity 
measurement lines (after the velocity component is sub- 
tracted of!) to the reference signal, the pressure shift can be 
found. Alternatively, one could use a single beam pointed 
directly across the flow (0=0, no velocity shift) as shown in 
FIG. 1 . This signal is illustrated in FIG. 5, where just the linear 
portion of the 1 f line shape is shown (dashed line on left side). 
From this signal and Eqn. (1), the pressure is 

P=P 0 +(A Vl Jb s ) (6) 

where Avp is the observed wavenumber shift between the 
reference pressure and flow signal . In a maimer similar to that 
above, the flow pressure in terms of measured voltages k and 
k' is 

P=P 0 MAvJ(2b,))*((k'+k)/(k'-k)) (7) 

As mentioned above, the mean shift of both crossed beams 
relative to the zero velocity reference signal could also be 
used to find the unidirectional pressure shift. 

For the purpose of illustrating the method of this invention, 
the slope of the If line shape near line center reflects the 
number density of the absorbing gas. By selecting an appro- 
priate absorption line tailored to an expected temperature 
range of interest, one can achieve a situation where the slope 
of the 1 f waveform near line center varies monotonically with 
density. The difficulty is that this density depends on both 
temperature and mole fraction of the species being measured. 


(fl-^)=Slope 1 y !< v D . 
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If one is working with oxygen in the inlet of a hypersonic 
engine, for example, the mole fraction of oxygen in air is 
always 0.2 1 , and the slope can be used to determine the local 
temperature directly. If one is measuring water vapor, there 
are two unknowns. However, there may be regions where the 
functional form of the slope versus density curve is such that 
the variation may be insensitive to either mixing ratio or 
temperature, so that one may be able to unravel both func- 
tions. One can easily make an additional measurement of the 
2f peak signals (just some minor additional electronics) that 
provides the peak absorbance and thus a second observable. 

In another example using the 1391.7 mn water vapor 
absorption line and assuming a pressure of 2 atmospheres and 
temperature of 1200 K, the fractional absorbance of combus- 
tor exhaust water vapor across an 8 cm path is approximately 

0.18. The shift in wavenumber with velocity for this line is 
approximately 0.002 cm -1 per 100 m/s. So for a Mach 1 flow, 
the shift would be about 0.008 cm - 1 . Pressure shifts for water 
vapor are -0.017 cm _1 /atm, so that a pressure of 2 atmo- 
spheres would have a shift (relative to a 0.25 atmosphere 
reference) of -0.030 cm -1 . 

Minimum detectable absorbance for systems using wave- 
length modulation spectroscopy is typically 10 -5 or better for 
a 1 Hz bandwidth. Since sensitivity scales as the square root 
of bandwidth, the expected sensitivity at 5 kHz would be 
7xl0~ 4 , equivalent to a SNR (signal to noise ratio) of 250. 
From these numbers, one predicts a measurement resolution 
of 0.8%. So for a velocity of 300 m/s, one could ideally 
resolve at 2.4 m/s. Of course, there is a wide variation in the 
possible mixing ratios, T, P, paths, and velocities that might 
exist in the engines being tested, but an expectation of better 
than 10 m/s fluctuation resolution at high bandwidth appears 
feasible. This precision should also apply to pressure, which 
at 2 atmospheres would be ±12 torr. 

In the preferred embodiment, and as readily understood by 
one of ordinary skill in the art, the apparatus according to the 
invention will include a general or specific purpose computer 
or distributed system programmed with computer software 
implementing the steps described above, which computer 
software may be in any appropriate computer language, 
including C++, FORTRAN, BASIC, Java, assembly lan- 
guage, microcode, distributed programming languages, etc. 
The apparatus may also include a plurality of such computers/ 
distributed systems (e.g., connected over the Internet and/or 
one or more intranets) in a variety of hardware implementa- 
tions. For example, data processing can be performed by an 
appropriately programmed microprocessor, computing 
cloud, Application Specific Integrated Circuit (ASIC), Field 
Programmable Gate Array (FPGA), or the like, in conjunc- 
tion with appropriate memory, network, and bus elements. 

Note that in the specification and claims, “about” or 
“approximately” means within twenty percent (20%) of the 
numerical amount cited. All computer software disclosed 
herein may be embodied on any non-transitory computer- 
readable medium (including combinations of mediums), 
including without limitation CD-ROMs, DVD-ROMs, hard 
drives (local or network storage device), USB keys, other 
removable drives, ROM, and firmware. 

Although the invention has been described in detail with 
particular reference to these preferred embodiments, other 
embodiments can achieve the same results. Variations and 
modifications of the present invention will be obvious to those 
skilled in the art and it is intended to cover in the appended 
claims all such modifications and equivalents. The entire 
disclosures of all references, applications, patents, and pub- 
lications cited above are hereby incorporated by reference. 


What is claimed is: 

1. A method for non-intrasively obtaining the local veloc- 
ity, pressure, and temperature of a gas, the method comprising 
the steps of: 

5 emitting two beams of laser light directed at two different 
angles across a flow of the gas; 
separately detecting the two beams; 
determining first harmonic line shapes from the detected 
beams using wavelength modulation spectroscopy and a 
to pair of single point measurements by a computer; 

calculating gas velocity and pressure from wavelength 
shifts of the line shapes by the computer; and 
calculating gas temperature from slopes of the line shapes 
by the computer; 

15 wherein wavelength of the two beams is alternated 
between two fixed wavelengths; and having an addi- 
tional signal channel to provide reference zero-velocity 
and reference pressure. 

2 . The method of claim 1 wherein a square wave is 
20 employed. 

3. The method of claim 2 wherein frequency of the square 
wave is at or greater than desired detection bandwidth. 

4 . The method of claim 1 wherein scanning full absorption 
line shapes is not required. 

25 5 . The method of claim 1 wherein gas velocity is deter- 

mined independent of gas density changes. 

6. The method of claim 1 wherein calculating gas tempera- 
ture comprises employing a mole fraction of the gas. 

7 . The method of claim 1 wherein calculating gas tempera- 
30 ture comprises measuring second harmonic peak signals. 

8. The method of claim 1 wherein the method operates at a 
bandwidth of up to about 40 kHz. 

9. The method of claim 1 wherein velocity is obtained at a 
resolution of about 10 m/s fluctuation or better. 

35 10. The method of claim 1 wherein pressure is obtained at 

about plus or minus 12 torr at 2 atmospheres. 

11 . .An apparatus for non-intrusively obtaining the local 
velocity, pressure, and temperature of a gas, said apparatus 
comprising: 

40 two emitters of beams of laser light directed at two differ- 
ent angles across a flow of the gas; 
detectors for the two beams; 

a computer determining first harmonic line shapes from the 
detected beams using wavelength modulation spectros- 
45 copy and a pair of single point measurements, calculat- 
ing gas velocity and pressure from wavelength shifts of 
the line shapes, and calculating gas temperature from 
slopes of the line shapes; 

a laser controller alternating wavelength of the two beams 
50 between two fixed wavelengths; and 

a reference cell providing reference zero-velocity and ref- 
erence pressure. 

12 . The apparatus of claim 11 wherein a square wave is 
employed by said laser controller. 

55 13 . The apparatus of claim 12 wherein frequency of said 

square wave is at or greater than desired detection bandwidth. 

14 . The apparatus of claim 11 wherein scanning full 
absorption line shapes is not required. 

15 . The apparatus of claim 11 wherein gas velocity is 
60 determined independent of gas density changes. 

16 . The apparatus of claim 11 wherein calculating gas 
temperature comprises employing a mole fraction of the gas. 

17 . The apparatus of claim 11 wherein calculating gas 
temperature comprises measuring second harmonic peak sig- 

65 nals. 

18 . The apparatus of claim 11 wherein said apparatus oper- 
ates at a bandwidth of up to about 40 kHz. 
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19. The apparatus of claim 11 wherein velocity is obtained 
at a resolution of about lOm/s fluctuation or better. 

20. The apparatus of claim 11 wherein pressure is obtained 
at about plus or minus 12 torr at 2 atmospheres. 



